The Supermassive Black Hole at the Heart of Centaurus A: 
Revealed by Gas- and Stellar Kinematics 

Nadine Neumayer^'® 

A European Southern Observatory, Karl-Schwarzschild-Str 2, 85748 Garching bei Miinchen, Germany 
B Email: nneumaye@eso.org 



Abstract: At less than 4 Mpc distance the radio galaxy NGC 5128 (Centaurus A) is the prime example 
to study the supermassive black hole and its influence on the environment in great detail. To model 
and understand the feeding and feedback mechanisms one needs an accurate determination of the mass 
of the supermassive black hole. The aim of this review is to give an overview of the recent studies 
that have been dedicated to measure the black hole mass in Centaurus A from both gas and stellar 
kinematics. It shows how the advancement in observing techniques and instrumentation drive the field 
of black hole mass measurements and concludes that adaptive optics assisted integral field spectroscopy 
is the key to identify the effects of the AGN on the surrounding ionised gas. Using data from SINFONI 
at the ESO Very Large Telescope, the best-fit black hole mass is M BH = 4.5ti;o x 10 7 M Q (from H 2 
kinematics) and Mbh = (5.5 ±3.0) x lO 7 M0(from stellar kinematics; both with 3a errors). This is one 
of the cleanest gas vs star comparison of a Mbh determination, and brings Centaurus A into agreement 
with the Mbh -a relation. 

Keywords: galaxies: individual (NGC 5128) - galaxies: kinematics and dynamics - galaxies: structure 
- galaxies: nuclei - techniques: spectroscopic 



1 Introduction 

Centaurus A (hereafter Cen A) is the posterchild of 
an active galaxy. Also known as NGC 5128 it is the 
closest elliptical galaxy, the closest recent merger, and 
hosts the closest active galactic nucleus (AGN). The 
distance to Cen A has long been und er debate and 
seems to converge to the value 3.8 Mpc (lReikuball2004l : 
iKarachentsev et al.| [2007: Har ris et alJ uOOQ). 

Cen A is one of the most studied galaxies in the sky 
and its proximity makes it an ideal candidate to test 
AGN models, the connection of merging and star for- 
mation, and the influence of the jet. Its proximity also 
makes it one of the prime targets to reliably measure 
the mass of the supermassive black hole, suspected 
to hide behind the obscuring dust lane. This promi- 
nent dust lane has hindered the study of CenA's black 
hole with optical HST spectroscopy as was common 
pract ice for other galaxies (e.g. M87. iMacchetto et al.l 
119971 1 and requires to move to the near-infrared wave- 
length for a detailed study of its properties. 
Interestin gly, the first measurements of Cen A' s black 
hole mass (|Marconi et al.ll200ll : ISilge et all2005h brought 
it al most a factor of ten above the so called M-a rela - 
tion (|Ferrarese fc Merritdl2000l : iGebhardt et alJl200Cf ). 
and made Cen A one of the largest outlier to this rela- 
tion. It was not clear, whether this offset to the M-a 
relation is intrinsic to Cen A, maybe due to the pres- 
ence of an AGN or its merger history, or whether it can 
be accounted for by the fact that seeing limited studies 
do not fully resolve the sphere of influence of the black 
hole, which would be at ~ 0'.'4, taking the predicted 
mass from the M-a relation llFerrarese fc Merrittll2000l : 
IGebhardt et al.ll2000l : iTremaine et al.l2002h forgranted . 



There are several methods to detect black holes 
in galaxy centres and to measure their masses. For 
nearby galaxies the most direct methods take gas and 
stars as kinematics traces of the nuclear region. Dy- 
namical models are then constructed to explain the 
kinematics. These models rely on a couple of assump- 
tions (circular gravitational motion of the gas, super- 
position of orbits for the stars). To test the reliability 
of both approaches it is important to compare their 
model results. In general, the gas kinematical method 
has the advantage of using a relatively simple mod- 
elling approach and dealing with relatively short ex- 
posure times. However, the gas kinematics might be 
influenced by non-gravitational motions that bias or 
even falsify the method. Moreover, not every galaxy 
nucleus has detectable gas emission lines. The stellar 
dynamical approach has the clear advantage that stars 
are always present in galaxy nuclei and that their mo- 
tion is always gravitational. However, the derivation 
of stellar kinematics requires long exposure times and 
gets very complicated in the close vicinity of an AGN, 
as the stellar absorption lines are diluted by the AGN 
continuum radiation. In addition, stellar dynamical 
models are v ery complex leading to pot e ntial indeter- 
mina cy (e.g. ICretton fc Emselleml [200i : IValluri et al.l 
2004) . Centaurus A is a prime example for testing both 
modelling approaches on our doorstep. It allows both 
gas and stellar kinematical studies at a spatial resolu- 
tion of ~ 2pc (corresponding to O'.'l) using either the 
HST or ground based adaptive optics instrumentation. 

This paper reviews the status of the measurement 
of the black hole mass in Centaurus A and shows how 
the advance in instrumentation is driving the field of 
black hole mass modelling. Section presents the dif- 
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ferent studies dedicated to the black hole mass mea- 
surement in Cen A, while Section [3] states what has 
been learned from these studies. Finally, Section [4] 
summarises the observational facts for Cen A's black 
hole. All studies on Cen A's black hole mass assumed 
a distance to Cen A of D=3.5 Mpc. The derived 
black hole mass directly relates to the adopted distance 
(Mbh oc D) and so all results need to be increased by 
8.5% given th e distance measurem ent of 3.8 Mpc rec- 
ommended bv lHarris et al.l l|2009t) . 

2 Dedicated Studies 

This section briefly lists and describes the studies that 
were performed to get the most accurate black hole 
mass determination to date of the supermassive black 
hole in Centaurus A. Each of these papers sets the 
state-of-the-art observations and modelling techniques 
available at the time the study was pursued. Gas and 
stellar kinematical studies are divided into two groups 
for comparative purposes. For details concerning the 
individual studies the reader is referred to the original 
papers. 

2.1 Gas kinematical studies 

Marconi et al. 2001 
Peering through the dust 

An extensive imaging study with WFPC2 and NIC- 
MOS on board the Hubble Space Telescope (HST) (PI: 
E. Schreier) has revealed the nucleus of Cen A in near- 
IR (K and H) and, for the first time, in th e optical (I 
and V ) at unprecedented spat ial resolution (Ma rconi et alj 
l200Cf ). ISchreier et~aH (|l998ft detected an elongated 
structure in Paa along the j et position ang l e and inter- 
preted it as a thin gas disc. iMarconi et alj (|200lh used 



first adaptive optics systems went online. Using NAOS- 
CON ICA (NACO) ilRousset et alj|2003t iLenzen et all 
Il998l ) at the ESO VLT. lHaring-Neumaver et all (1200(1) 
(herea fter HN+06) followed the approach of lMarconi et alJ 
(|200ll ) and used long-slit spectra at the same position 
angles to probe the kinematics of the ionised gas. Due 
to the correction of the atmospheric distortion, the 
spatial resolution of the spectra increas ed by almost a 
factor of five compared to the study of lMarconi et al.l 
(|200lh and provided images and spectra at (or almost 
at) the diffraction limit of the VLT (0'.'06). The data 
points are shown in Figure Q] along with t he [Fell] kine- 
matics derived by IMarconi et alj (120011 ) from ISAAC 
seeing limited data. HN+06 note that the [Fell] veloc- 
ity dispersion is very high along all slit positions and 
cannot be fully accounted for by unresolved rotation. 
They have thus included a pressure term in their gas 
kinematical model. Moreover, they take into account 
the underlying potential of the stars (that was not ac- 
counted for by Marconi et al.l (|200ll )) and model the 
[Fell] kinematics in the joint potential of the stars and 
an additional point mass, assumed to be the black hole. 
Their data along four long-slit position angles does also 
not constrain the inclination angle of the modelled gas 
disc well. Their best fitting model gives i — 45° and 
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The decrease in black hole 
is due to (i) 



mass compared to IMarconi et al.l (|200ll 
accounting for the stellar potential, (ii) the higher reso 
lution of the data, and (iii) the higher inclination angle 
of the gas disc. 



the n ear-infrared spectrograph and imager ISAAC (|Moorwood e t at,' 
19991 ) on the ESO Very Large Telescope to retrieve 
spectra in J, H, and K s Band under excellent seeing 
conditions of 0'.'5 - 0'.'6. They used the HST narrow 
band images as a reference to position the slit along 
three position angles and derived the kinematics for 
both ionised and molecular gas from their emission 
lines (Pa/3, [Fell], Br7, and H2). The ionised gas kine- 
matics are modelled by a thin disc model in the grav- 
itational potential of a point mass, representing the 
black hole. Free parameters in the model are the incli- 
nation angle of the modelled disc, the black hole mass 
and the systemic velocity of the system. The inclina- 
tion angle and the central mass are strongly coupled, 
since the amplitude of the rotation curve is propor- 
tional to \/Mbh x sin(i). The ISAAC data do reject 
disc inclinations > 60°, but do not constrain the disc 
inclination further. This gets reflected in the large un- 




Figure 1: Velocity curves extracted from the 
[Fell] emission line at A1.644^m from NACO (di- 
amonds) and ISAAC data (triangles), along three 
slit positions centred on the nucleus of Cen A 
(P.A.= -44.5°, 32.5°, and 82.5° from left to right). 
The NACO data have almost factor of five better 
resolution than the ISAAC data and are overall 
consistent. 



certainties of the black home mass determination. The 



best fit value is Mbh 



1+3.0 



x 10 8 M Q 



Marconi et al. 
Pushing STIS 



2006 

to the reddest end 



Haring-Neumayer et al. 2006 
Adaptive Optics works! 

In the year 2004 the resolving power of ground based 
instrumentation has taken a big leap ahead when the 



Knowing exactly where to look, IMarconi et all ( 2006) 
have pointed the HST spectrograph STIS on the dust 
enshrouded nucleus, that is not visible in the optical 
wavelengths. Pushing STIS to the reddest limit, they 
used the [SIII] A9533 A line to study the kinematics of 
the ionized gas in the nuclear region with a 0'.' 1 spatial 
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resolution. The STIS data were analysed in conjunc- 
tion with the gr ound-based ne a r-infr ared VLT ISAAC 
spectra used bv lMarconi et all l|200lf ). and have a gain 
in spatial resolution of almost a factor of five (see Fig- 
ure [2]) • The gas kinematical analysis provides a mass 
of Mbh = (1.1 + 0.1) x 10 8 M Q for an assumed disc in- 
clination of i = 25° and Mbh = (6.5±0.7) x 1O 7 M for 
i = 35°, the larges t inclination va l ue all owed by their 
data. In this study, [Marconi ct al. (2006) did extensive 
tests on the influence of the choice of the parametri- 
sation of the surface brightness profile of the assumed 
gas disc, that is used to luminosity-weight the mod- 
elled emission lines and to resemble the observations. 
They found that the associated systematic errors are 
no larger than 0.08 in logMsH comparable with statis- 
tical errors and indicating that the method is robust. 
However, the intrinsic surface brightness distribution 
has a large impact on the value of the gas velocity dis- 
persion. They conclude that a mismatch between the 
observed and modelled velocity dispersion is not nec- 
essarily an indication of non-circular motions or kine- 
matically hot gas, but is as easily due to an inaccu- 
rate computation arising from too course a model grid, 
or the adoption of an intrinsic brightness distribution 
which is too smooth. Therefore, they did not include 
pressure support in their model. The stellar potential 
is accounted for by the model. The spatial resolution 
of the data is ~ O'.'l the same as in HN+06 and both 
studies use ionised gas as the kinematic tracer ( [SIII] 
and [Fell], respectively). The largest discrepancy be- 
tween the two studies is the modelled disc inclination 
angle. iMarconi et al.l l|2006h find that the inclination 
angle is i < 35°, while HN+06 find their best model 
at i = 45° and exclude smaller inclination angles on 
the basis of the jet inclination angle (assuming a or- 
thogonal jet-disc geometry and a jet inclination ang le 
of~50° |Tingav et aill 19981 ; lHardcastle et al]|200a i). 
The inclination angle remains the stro ngest source of 
uncert ainty in the model. The results of lMarconi et al.l 
(2006) and HN +06 are in very good agreement, and 
also agree with Marconi et al.l l|200l] ) inside the error- 
bars. 



Krajnovic, Sharp & Thatte 2007 
Integral Field Spectroscopy to the rescue! 

All previous approaches were restricte d to a limited 
numb er of slit position angles (3 for | Marconi et al.l 
(|200ll 1. 4 for HN+06, and 2 for lMarconi et all (|2006l V) 
and did not allow to trace the kinemat ics of the gas 
outsid e these stripes. This changed when lKrainovic et al] 
(2007) analysed the data from the Cambrid ge Infrared 
Pano ramic Survey Spectrograph (CIRPASS. lParrv et al] 
2004), an integral field spectrograph mounted on the 
Gemini South telescope. They detect two dimensional 
distributions of the following emission lines: [PII], [Fell], 
and Pa/3, and extracted spatially resolved 2D kinemat- 
ics from Pa/3 and [Fell]. All emission- line regions are 
part of the same kinematic structure which shows a 
twist in the zero-velocity curve beyond ~ 1" (for both 
Pa/3 and [Fe II]). The kinematics of the two emis- 
sion lines are similar, but the Pa/3 velocity gradient 
is steeper in the centre while the velocity dispersion is 
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Figure 2: Comparison of average velocities from 
STIS derived from [SIII] emission and ISAAC de- 
rived from [Fell] and Paa, showing the increase in 
spatial resolution of almost a facto r of five. This 
i s a r eproduction of Fig. 8 from IMarconi et al.l 
(2006) 



low everywhere. The velocity dispersion of the [Fell] 
emission is relatively high featuring a plateau, approx- 
imately oriented in t he same way as t he cen tral part 
of the warped disc. IKrainovic et al.l (|2007f ) use the 
2D kinematic information to test the hypothesis that 
the ionized gas is distributed in a circularly rotating 
disc. Assuming simple disc geometry they estimate the 
mass of the central black hole using Pa/3 kinematics, 
which is consistent with being distributed in a circu- 
larly rotating disc. Their best fit model gives Mbh = 
8.2514 25 x 1O 7 M , for PA = -3° and i = 25°, and a 
stellar mass-to-light ratio of 0.0 M/Lq, indicating that 
the Pa/3 gas disc is fully under the gravitational influ- 
ence of the black hole. To test their model against pre - 
vious gas dynamical models, IKrainovic et al.l (|2007t ) 
constructed models with different geometrical orienta- 
tions taken from previous studies. They find that these 
give similar results, and are often statistically indis- 
tinguishabl e. However, the mass predicted by the M- 
g rel ation ijFerrarese fc Merritt 2000; Gcb hardt et al] 
2000) is excluded by all their models. 

Neumayer et al. 2007 

Integral Field Spectroscopy combined with 
Adaptive Optics 

Another advance in observing techinque allowed a big 
jump ahead in the studies of black hole mass mod- 
elling. Combining the resolving capabilities of adap- 
tive optics with the power of int egral field spectroscopy , 
makes SINFONI o n the VLT (|Eisenhauer et al] 12003 : 
iBonnet ~et al. 2004) the perfect instrument for this case. 
SINFONI works in the near-infrared and allows to peer 
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through the dust lane hiding Cen A's nucleus. The 
adapt ive optics assisted data presented bv lNeumaver et all 
(2007) have a spatial resolution of 0'.'12 and cover the 
central 3" x 3". The data cubes in H- and K-band 
exhibit a wealth of ionised gas emission lines: [Fell], 
[SiVI], Hel, Br7, [CaVIII], and several transitions of 
H2. In addition the data enable the extraction of stel- 
lar kine matics (which will be pre sented in detail in Sec- 
tion [2]2];[CaE2enirreFal][200i). The high spatial res- 
olution two-dimensional maps of the ionised gas show 
clearly two distinct kinematic features: rotation and 
transverse motion (see Figure [3} . The transverse mo- 
tion is oriented along the jet direction (PA=55°) and 
is strongest for the high-ionisation lines ([SiVI] and 
[CaVIII]). However, also the tracers used in previous 
studies like [Fell] and Br7 are influenced by this non- 
gravitational motion that has not been accounted for 
by the dynamical models. Long-slit spectroscopy was 
not able to spot this additional kinematic feature and 
therefore leads to an over-simplified picture, given the 
assumption that the gas moves in a circularly rotating 
disc, which underlies the gas dynamical models. 
The picture is surprisingly different when looking at 
the kinematics of molecular hydrogen rather than the 
ionised gas. This kinematic tracer shows no kinematic 
component along the jet direction and seems to be 
dominated by gravita tional motion only. Therefore, 
iNeumaver et al.l l|2007h have based their dynamical model 
on these kinematics. The two dimensional kinematic 
maps constrain the geometry of the disc much better 
than the few long-slit position angles have done be- 
fore. The disc is not flat, but appears warped. The 
disc parameters, i.e. the inclination angle and the ma- 
jor axis position angle, were fitted usin g the kineme- 
try software by lKrainovi(5 et alJ l|2006h . The fitted 
mean inclination angle of the H2 gas disc is 45° ± 
7° and the mean position angle of the major axis is 
155°. The inclination an gle is in very good agree ment 
with the value derived by Hardcas tle et al.l (|2003f ) from 
VLA data (20° < i < 50°), it is, however, somewhat 
smaller than the v alue from VLBI data derived by 
iTingav et all (|l998fl (50° < i < 80°). This well con- 
strained disc parameters enable a significant reduction 
of the uncertainty of the modelled black hole mass. 
The assumed mass-to-light ratio of the underlying stel- 
lar body, M/L K = (0-72 ± 0.04) M R /L R , is taken from 
the stellar ki nematical model bvlSilge et al.l ll2005h and 
confir med bv lCappellari et alJ (|2009h . INeumaver et al.1 



CIRPASS/Gemin 



(2007) modelled the H2 kinematics on tilted rings fol- 
lowing the kinemetry model, keeping the relative incli- 
nation angle between the rings fixed but varying the 
overall inclination. The overall best-fit black hole mass 
is Mbh = 4.5l*'o x 10 7 M Q and the best fit mean incli- 
nation angle is 34° ±4° (3 a error bars) . 



2.2 Stellar kinematical studies 

Silge et al. 2005 
Near-Infrared long-slit spectra 

Using the Gemini Near-Infrared Spe ctrograph (GNIRS) 
on Gemini South. [Silge et all l|2005h measured the stel- 
lar kinematics under very good natural seeing condi- 
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Figure 3: The [Fell] velocity maps compared 
from SINFONI adaptive optics and CIRPASS 
seein g limited data. Th e Fi gures are taken 
from I Neumaver et al.l ( 2007[ ) and iKrainovic et al.l 



(|2007r ) respectively. The overall kinematic struc- 
ture of the two maps is in good agreement. How- 
ever, the adaptive optics assisted SINFONI map 
shows detailed kinematic substructures and reveals 
the influence of the jet on the kinematics. 



tions of 0'.'45 and 0'.'6. They put the slit along two posi- 
tion angles across the nucleus of Cen A, centred on the 
AGN. The first observation was made with the slit ori- 
ented perpendicular to the inner dust ring (along the 
major axis at large radii), while the second observation 
was made with the slit oriented parallel to the dust 
lane but offset from the centre by / /85. The stellar 
kinematics are extracted using the region around the 
CO band heads at 2.3 /im. The central spatial bins of 
the first observations are dominated by emission from 
the AGN and the CO band heads are diluted by this 
emission. To recover the stel lar signatures inside the 
central 0'.'3. ISilge et all (120051 ) have removed the AGN 
emission by measuring the AGN contribution by its di- 
lution of the equivalent width of the (2— l) 12 CO band 
head, assuming that the equivalent width is largely 
constant with radius in this galaxy for regions out- 
side the AGN. After subtracting the AGN contribu- 
tion, the stellar kinematics were extracted from the 
spectra by fitting a stellar template spectrum (con- 
volved with a velocity profile) to the observations. The 
choice of t emplate star is importan t in order to avoid 
mismatch (|Silge fc GebhardtJ 120031 ) and so a variety 
of template stellar spectra is simultaneously fitted for 
the velocity profile and the stellar template weights. 
As a result the fitting procedure provides the line-of- 
sight velocity distribution (in a non-parametric form) 
as well as the stellar population information. The lu- 
minosity weighted stellar velocity dispersion along the 
slit parallel to the dust disc is 138 ± 10 km s _1 . 
The black hole mass is determined using axisymmet- 
ric orbit-based models. The full line-of-sight veloc- 
ity distribution is fitted to the data and the best fit 
parameters, stellar M/L and Mbh , are derived via a 
y 2 minimisation (the interested reader is referred to 
ISilge et all (|2005l 'l for details) . From the long slit data 
along two slit positions alone it is not possible to iden- 
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tify the rotation axis and so the modelling procedure 
was performed twice: first matching the rotation axis 
to the dusk disc and second matching the rotation axis 
to the photometric major axis and non-rotating kine- 
matics along the dust disc. It turns out that the two 
assumptions give very similar black hole masses but 
the second approach results in significantly lower \ 2 
values. 

Although NGC 5128 appears very round, its actual 
shape is unknown a nd it could be int rinsically quite 
flattened. Therefore, ISilee et all (|2005h repeated their 
models for three different intrinsic shapes (i.e intrinsic 
axis ratios) of Cen A: intrinsically round (1:1:1), at an 
inclination of i = 20° (1:1:0.5), and i = 45° (1:1:0.9). 
They find that this drastic change in flattening does 
not change the black hole mass much. The \ 2 of the 
intrinsically round model is significantly less than ei- 
ther flattened model, and so this is adopted as the 
best-fit model with M B h = 2.4lo 2 x 10 8 M Q and stel- 
lar M/L K = 0.1 
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Cappellari et al. 2009 

Integral field spectroscopy combined with 
adaptive optics 

Using the adaptive optics assisted integral field spec- 
trograph SINFONI at the ESO Very Large Telescope, 
Cappellar i et alj (|2009l ) extracted the stellar kinemat- 
ics from the nuclear region of NGC 5128 with a spatial 
resolution of 0'.'17. This means an increase in spatial 
resolution by a facto r of three to four c ompared to the 
seeing limited data of lSike et ail (|2005T) . The high S/N 
data cover the central 8" x 8" and allow the extraction 
of the stellar line-of-sight velocity distribution, which 
is given in parametric form of a Gauss-Hermite expan- 
sion. At radii where the nucleus dominates, the stellar 
absorptions are diluted, resulting in a strong decrease 
in the observed line-strength 7. An inaccurate mod- 
elling of the non-thermal dilution can cause large errors 
in the measured velocity dispersion er. Similarly, the 
rise in the non-thermal continuum could be incorrectly 
interpreted as a variation of the stellar population, re- 
quiring a change in the stellar template mix and also 
producing an error in a. As there is no evidence for 
a sudden change in the population in the nucleus of 
Cen A, the authors argue the safest choice is to as- 
sume that the stellar template is fixed and to model 
the non-thermal continuum via additive polynomials. 
This approach allows the a to be reliably extracted in 
the high-resolution observations down to > 0'.'2. For 
the first time, a low-level clear sense of stellar rotation 
is detected for Cen A which is counter-rotating with re- 
spect to the gas. The global kinematical major axis is 
at PAkin = 167 ± 8° which is in between the dust disc 
and the galaxy's photometric major axis, and there- 
fore b etween the two position angles that Silgc ct al. 
(2005) assumed for their kinematic models. The de- 
rived systemic velocity is v Byst = 531 ± 5 km s -1 . 

Having in mind the complex dynami cal and mor- 
pholog ical structure that Cen A possesses. ICappellari et ail 
(2009) constructed two types of models to test the sen- 
sitivity of the Mbh estimate to the assumed geometry 



Figure 4: Spherical anisotropic Jeans model to 
check the Mbh determinations. The root mean 
square stellar velocity, v rms = \/v 2 + er 2 , is plotted 
as a function of radius from the centre of Cen A, 
derived from SINFONI observations (100 mas, red 
open circles; 250 mas scale, b lue open squares) 
and from GNIRS observations of lSike et all (|2005h 
along the major axis (filled downward triangle) 
and the minor axis (filled upward triangle). The 
discrepancy between the SINFONI and GNIRS 
data inside 1" is due to a different subtraction 
of the AGN continuum. The solid line shows 
the prediction of an anisotropic Jeans model, hav- 
ing the best-fitting Mbh and M/L derived from 
the Schwarzschild model, the dashed lines rep- 
resent the upper and lower 3tr errors in the 
Schwarzschild model. The dotted lines has Mbh = 
2 x 10 8 M ff , co r respo nding to the best-fitting value 



of ISilge et afl (l2005lh This figure is reproduced 



from ICappellari et al.1 (2009) 



and dynamics; (i) an axisymmetric orbit-based model, 
to reproduce in detail the kinematic observations and 
(ii) a simple anisotropic spherical Jeans model to qual- 
itatively check the Mbh determination. This test gave 
very consistent results and is shown in Figure U along 
with a comparison to the data presented bv lSilge et al.l 
(2005). The best-fitting value for the black hole mass 
from the Schwarzschild model is Mbh = (5.5 ± 3.0) x 
lO 7 M (3o- errors) and the corresponding best fitting 
stellar M/L K is (0.65 ± 0.15)M/L Q . 



3 Lessons learned 

The ever increasing accuracy of the black hole mass 
measurement in Centaurus A over the past eight years 
reflects the advancement in observing techniques, es- 
pecially in the near infrared. High spatial resolution 
observations are crucial to determine the influence of 
the central black hole on the gas and stellar kinematics. 
Even in an object that is as nearby as Cen A, where 
the sphere of influence of the black hole is comfortably 
resolved (Rbh ~ OVj, using the mass measurement of 
ICappellari ela l. (2009)) using HST or adaptive optics 
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observations. 

Adaptive optics actually works! With the advent of 
near-infrared adaptive optics assisted spectrographs, 
such as NACO at the VLT, even dust enshrouded galaxy 
nuclei became accessible at a spatial resolution of O'.'l. 
This gets even more powerful when combined with 
integral-field spectroscopy. Mapping the gas and the 
stars in 3D allows to directly and simultaneously com- 
pare the morphology and kinematics of different gas 
species plus the stars. Having this powerful tool in 
hand, the influence of the inner jet on the kinemat- 
ics of the ionised gas in Cen A could be revealed, and 
moreover, molecular hydrogen could be identified as 
the ideal gas tracer for the central gravitational po- 
tential. The physical state of the gas is therefore very 
important when using it as a tracer for the dynam- 
ical models. The de cr ease of the bes t -fit v alue from 
iMarconi et"al] (|200lT ) to lMarconi et alj (|2006T ) and HN+06 
was due to the increase in spatial resolution plus taking 
into account the stellar contribut ion to the po t ential , 
while the decr ease in Mbh fromlMarconi et al.l (J2006) 
and HN+06 to iNeumaver et alJ (120071 1 was due to the 
fact that the kinematic tracer changed from ionised 
gas to molecular hydrogen. 

The presence of an AGN can definitely influence the 
kinematics of the gas, while the stellar kinematics should 
be unchanged by this. However, the extraction of the 
stellar kinematics from the spectral absorption fea- 
tures gets increasingly difficult in the close vicinity of 
the AGN, as the non-thermal continuum dilutes the 
stellar absorption lines and needs to be accounted for 
in the extracti on process. This is the main difference in 
the an alysis of lSilge et alJ 1120051) a nd lCappellari et afl 
<|2009T ). While ISilge et al.l $2063 ) first subtract the 
AGN contribution an d then fit the stell a r line -of-sight- 
velocity distribution, ICappellari et all (|2009ft include 
the fit of the AGN continuum in the the extraction 
of the stellar kinematics. To get reliable kinematics 
high signal-to-noise data are required. This is a very 
interesting lesson that we learned from Cen A, and we 
should be cautious when extracting kinematics from 
other, more distant objects. Cen A is indeed the clos- 
est AGN and at the same time it is very complex. 
Every leap in instrumentation development is likely to 
reveal more complex substructures. This warrants our 
continuous attention, in order to reveal intrinsic prop- 
erties in the data and understand shortcomings in the 
models that aim to predict the observations. 



4 Summary 

In principle all presented Mbh determinations for Cen A 
using gas kinematical models agree within the error 
bars. The only significant disagreement in Mbh is with 
the previous stel lar dynamical Mbh determination by 
Silge et al.l ll2005[) and subsequent measurements (HN+06, 
Marconi et al 2006: lKrainovic et al l2007l : lNeumaver et alJ 
2007l ; ICappellari et all2009ft . The disagreement is likely 
caused by a difference in the data quality and in the 
treatment of the contribution from the central non- 
thermal continuum in the kinematic extraction (com- 
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pared to ICappellari et al]|2009h . It is not due to dif- 
ferences in the details of the modelling methods. 
The latest mass determination for the black hole at 
the nucleus of Centaurus A, using high resolution inte- 
gral fiel d stellar kinematics, g ives Mbh = (5.5 + 3.0) x 
10 7 M^ (|Cappell ari et al . 200§) and agrees very well with 

the value Mbh = x 10 7 M B determined using H2 

gas k inematics from the same SINFONI data ( Ncumavcr et al. 
2007) . This gas versus stars Mbh com-parison consti- 
tutes one of the most robust and accurate ones, due 
to a very well resolved black hole sphere of influence 
{Rbh ~ 0'.'7 compared to a PSF FWHM~ 0'.'12) and 
thanks to the use of high-resolution integral-field data 
for both kinematical tracers. 

Compar able and equa l ly suc cessf ul comparisons wer e 
done bv lShapiro et al.l (|2006h and lSiopis etHI (I2009D. 
while a less good agreement was found in lCappellari et al] 
(2002), likely due to the disturbed gas kinematics. There 
are certainly weaknesses in both the stellar and gas 
dynamical models, as e.g. the assumption of circular 
gravitational motions for gas kinematics and the com- 
plexity for the modelling of stellar kinematics. How- 
ever, the agreement between the two modelling ap- 
proaches for a growing number of cases strengthens 
the reliability of both methods. 

With the latest gas and stars Mbh determination 
Cen A lies wit hin the errors on the M b h — cr relation as 
given by ei therlTremaine et al.l (l2002h . lFerrarese fc Ford! 
(2005), or lGiiltekin et al.l l|2009h . and is in agreement 
with the re lation of black hole mass to bulge mass 
as given by lHaring fc Rixl l|2004h (see Figure • As 
the value of Cen A's black hole mass has changed 
since 2001, so has the Mbh — a relation. The scat- 
ter in the relation seems to depend on the Hubble 
type of the c hosen galaxies and is smallest for ellip- 
tical galaxies (|Hull2008l ; rGrahamll2008l ; lGultekin et all 
120091 ). The latest SINFONI M B h measurement is in 
agreement with the relation for ellipticals (Figure 5). 
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Figure 5: Cen A's black hole mass measurements are plotted on top of t he black hole m a ss gal axy scaling 
relations. The left panel shows the Mbh -o relation reproduced after iGiiltekin et al with the 

overplotted line being the best fit relation to the elliptical galaxy subsample of all plotted galaxies. The 
right panel shows the Mbh -Meuigo relation as presented in lHaring fc Rbd ( 2004) with Cen A values over- 
plotted. The Cen A plot symbols as well as their time sequence are indicated in the upper left corner of 
the left and right panels, respectively. Values published in the same year are separated horizontally for 
clarity. In the left panel (and for Cen A measurements) triangles refer to gas kinematical measurements, 
while circles refer to dynamical models using stellar kinematics. 



